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SUMMARY 


Williams type modal solutions of the elementary and Timoshenko beam 
equations are presented for the transient response of several uniform 
beams to a general applied load. Example computations are shown for a 
free-free beam subject to various concentrated loads at its center. The 
discussion includes factors influencing the convergence of modal solu- 
tions and factors to be considered in a choice of beam theory. Results 
obtained by two numerical procedures, the traveling-wave method and 
Houbolt's method, are also presented and discussed. 


INTRODUCTION 


The problem of obtaining the response of elastic structures to 
rapidly applied loading is of continuing concern to the aircraft indus- 
try inasmuch as aircraft structures must withstand blasts, landing 
impacts, and a variety of other transient loads. In order to study the 
various factors involved in this problem. It is desirable to consider 
simplified structures for which thorough studies are possible. Among 
the simplest examples of continuous elastic structures are uniform beams. 
Consequently, beams have been the subject of a considerable number of 
transient response investigations, and a variety of solutions of partic- 
ular beam problems are scattered throughout the existing literature. 

(See, for example, refs. 1 to 7» For an extensive bibliography, see 
ref. 7-) 


It Is the purpose of the present paper to provide a relatively com- 
plete source of useful modal solutions and to discuss the factors influ- 
encing the convergence of modal solutions and factors involved in the 
choice of the proper beam theory to be used In an analysis. To this, end, 
a consistent presentation is made of Williams type modal solutions for 
the response to a completely general transient load of three pertinent 
uniform beams (a free-free beam with a concentrated mass as its center, 
a cantilever beam, and a simply supported beam). (Some duplication of 
the existing literature is included for completeness.) Solutions, based 
on both the elementary and Timoshenko beam theories, are obtained by a 
process which can be readily extended to the solution of problems with 
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time -dependent boundary conditions. The application of the method is 
illustrated for the case of the free -free beam with a concentrated mass, 
and results for a! 1 the beams are summarized in tables I and II. In 
addition, some typical computed results are shown for a free-free beam 
subjected to various concentrated loadings. 

Another purpose of the present paper is a critical discussion of 
two numerical procedures, the traveling -wave method (ref. 6) and Houbolt' 
method (ref. 8). Ihe procedures are briefly described and computations 
made with both methods are compared with the modal results. 


SYMBOLS 


"s 

C 

C 1 


E 

f 

f 

G 

I 

k 

^RI 


effective shear-carrying area of cross section 
arbitrary -constant 

propagation velocity of bending discontinuities, 

I mr 2 


propagation velocity of shear discontinuities. 
Young' s modiolus of elasticity 
applied concentrated load 


dimensionless applied concentrated load, ^ 

shear modulus of elasticity 
moment of inertia of crosB section 
integers 


fl 2 


2,/m 


dimensionless frequency parameter, ct )Z 


dimensionless rotary inertia parameter, r/Z 


V m 


dimensionless transverse shear parameter. 




length of beam (half-length in case of free-free beam) 


Z 
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M 

M 

M s 

“i 

m 

me 


■bending moment (see fig. 1) 
dimensionless 'bending moment, Ml/EI 
dimensionless static bending moment 

generalized mass 
mass per unit length 
concentrated mass 

ratio of the concentrated mass to total mass of the beam, 
m c /ml 

generalized force 

applied distributed load (see fig. l) 
dimensionless applied distributed load, 
cross-sectional radius of gyration 
time 

transverse shear force (see fig. 1) 
dimensionless transverse shear force, VI^/eI 
dimensionless static transverse shear force 

coordinate along the beam 
deflection (see fig. 1) 
dimensionless deflection, y/l 

dimensionless translational component of ith natural mode 
dimensionless static deflection 


y r dimensionless rigid-body translation 
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Pi “ * 1 ^ ( k s 2 + ^I 2 ) + j|( k s 2 “ %I 2 ) + ~2 


r I = „ 2 


a ± 2 + k ± 2 k s 2 


Pi" - ki^s 2 


8 ( 1 ) 


/ r“ 

Dirac delta function l&(f-) = 0 for | ^ 0; / &({-)d| = 




0i 

* 

*i 

t c 


dummy variable of integration 

Pi 

sin - _±. sinh 

cos + y- cosh 

dimensionless space coordinate, x/l 
dimensionless time, - , 

1 2lm 

ith generalized coordinate 
rotation of beam cross section 
rotational component of ith natural mode 

static rotation of beam cross section 


co 

H r ) 


circular frequency of natural vibration 
step function (j( t) = 0 for t < 0j J( t) 
Matrix notation: - 
[ J rectangular matrix 

L J row matrix 

| | column matrix 

£ J diagonal matrix 


= 1 for t > 




t 
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Primes and Roman n ume ral superscripts are used to denote partial dif- 
ferentiation with respect to § . Dots denote partial differentiation 
with respect to r. 

WILLIAMS TYPE MODAL SOLUTIONS 


In normal-mode solutions for the response of beams to transient 
loads , the response Is expanded in terms of a series of normal modes of 
the beam. The coefficients of the expansion (the generalized coordinates) 
are determined from the governing differential equations and the boundary 
and initial conditions. Williams type modal solutions (ref. 2) differ 
from, ordinary normal-mode solutions by virtue of the Isolation of that 
portion of the response which may be obtained in closed form by a process 
of direct integration - the so-called "static" portion of the response. 
Only the remaining "dynamic" portion of the response is expanded in 
series form. 

The advantage of the Williams method over ordinary modal solutions 
is its ability to yield, for many loading conditions, a more accurate 
result with the same number of terms in the series. (See, for example, 
refs . 4 and 5 • ) It is particularly advantageous where the response func- 
tion is discontinuous. (An example of this is the determination of the 
shear due to a concentrated load.) The discontinuity is contained 
exactly in the separated static portion of the response and the series 
is only required to reproduce a continuous remainder. 

In the Williams method, the isolated portion of the response is 
termed static because significant parts of the inertia forces are ignored 
in its determination. In general, however, it is time dependent by virtue 
of the time dependence of the applied load and of the nonhcmogeneous time- 
dependent boundary conditions if such are imposed. In the case of beams 
with a fixed point of reference, such as cantilever or simply supported 
beams, all Inertia forces are ignored In the determination of this static 
part of the response j for beams with rigid-body freedoms, however, the 
inertia forces due to the rigid -body motion must be taken into account. 

One method of obtaining Williams type modal solutions is illustrated 
herein for both the elementary and Timoshenko beam theories. 


Elementary Beam Theory 

Basic equations .- The motion of a beam subjected to an applied load 
of intensity q(x,t) Is usually taken to be governed by the Bernoulli - 
Euler equation 
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b 2 y b 2 y 

—4 + m = 


EI^ + 
dx 2 


dt £ 


<1 


( 1 ) 


where x is the cpordinate along the beam, t is time, y(x,t) is the 
deflection (see fig. 1), El(x) is the bending stiffness of the beam, 
and m(x) is its mass per unit length. The internal bending moment M(x,t) 
and the shear force V(x,t) at any cross section (see fig. 1) are given by 


M 



and 


V = 




EI 


dfy 

bx.2- 


( 2 ) 


(3) 


respectively. 

For a uniform beam, these equations may be written in the dimension- 
less forms 


where 


y 17 + y - 5 


M = -y" 



y(|;T) 

S(|,t) 

V(|,t) 


_ y 

i 

= 

EI 

. 

' EI 


5(6, t) -g! 


7 

(5) 

(6) 


*■ 


f 
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and l is the length of the beam, or half-length in the case of a free 
free beam. The primes denote partial differentiation with respect to 


l 


— and dots denote partial differentiation with respect to r 


l2' m 


Symmetrical free-free beam with concentrated mass .- For symmetrical 
motion of a uniform beam having free ends at | = 1 and | = -1, atten- 
tion is restricted to the portion 0 Is | ^ 1 with boundary conditions 
stated in the form 


y' (0 ,t) = 0 


(7a) 


y'"(0,T) = 0 (7b) 

y"(l,T) = 0 (7c) 

y"'(l,T)=0 (74) 

If, in addition, the free-free beam has a concentrated mass 2m c located 
at the center | = 0, the influence of this mass may be introduced into 
the problem by changing the boundary condition, equation (7b), to 
y ,n (0 ,t) + m c y(0,T) = 0 where m c = m c /mZ . On the other hand, the bound- 
ary conditions, equations (7), may be left unchanged and the differential 
equation (eq. (4)) altered to 


y 17 + [l + m c 6(|)] y - 5 (8) 

where 5(|) is the Dirac delta function. In the solution that follows, 
the latter alternative is chosen. 

The b eam is assumed to be initially at rest and undeflectedj that is. 


y(S,o) = y(!,o) = 0 (9) 

Then the response to a general symmetrical load q(i,r) may be obtained 
in the Willi ams form by the following procedure. 


Hie solution is assumed in the form 
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y(l,*r) = y r (r) + y e (|,T) 0 ± (t) y ± (4) 

i=0 


9 


( 10 ) 


The quantity y r is the rigid-body translation of the free-free beam. 
It is determined to satisfy the differential equation 

y r ( T ) - ■ 1 ' , f a(i,T)d| (11) 

1 + Be Jo 

and the initial conditions 


y r (o) = y r (o) = 0 (12) 

The quantity y s (E, ,t) is the static deflection determined to satisfy 

= q(|,r) - [l + m c &U)]y r (T) (13) 

and the cantilever boundary conditions 


y s (0,T) =0 
y s ’ (o,t) = 0 
y s "(i,T) = 0 

y m (i,T) = 0 


(ik) 


Note that, by virtue of the definition of y r , y g also satisfies the 
boundary condition y B '*'(0,T) = 0. Finally, the shapes y ± ) (where 
i = 0, 1, 2, . . .) are the natural vibration modes of the beam satisfying 


y. (4) 


1 + m c 5(4)] hi 2 3^(4) 


9 

# 


( 15 ) 



7 


MCA TN b2kb 


9 


and the boundary conditions 


ft' ( 0 ) = 0 " 
yi m (o) = o 
yi'(i) - o 
yi ,M CD = o 


(16) 


<» 


* 


where the dimensionless frequency coefficients 

i = 0, 1, 2, . 


. ) are defined by kj_ = cjo^I y — . 

V El 


(where 

Further, it can be 


shown that the modes y^(0 satisfy the orthogonality condition 

,1 

0 


J 1 [l + I c 8(g)]y ± (* ) yj (| )d| = 0 (i / i) 


(IT) 


Note that, by virtue of the arbitrary selection of a datum plane for y 8 , 

the dynamic portion of the response, in general, still contains a rigid- 
body component (i = 0) . As defined, the total deflection y(g ,t) satis- 
fies the boundary conditions (eqs . (7))* There remains the problem of 
determining the coefficients 0^(t) so that the differential equation 

(eq. (8)) and initial conditions (eqs. ( 9 )), are satisfied. 

If expression (10) is substituted into differential equation (8) and 
equations ( 15 ) and ( 15 ) are taken into account, the differential equation 
is reduced to 


fcoo + k^iC-oJi^Ci) = -y B (|,r) (18) 

i=0 

Multiplying equation (18) by [jL + £<, 6(|)Jyj(|) and integrating with 

respect to | from 0 to 1 yields, in view of equation ( 17 ), the 
following result: 


I 


* 


$j_( T ) + (t) 


P±(t) 

mi 


\ 


(i = 0, 1, 2, . . .) ( 19 ) 
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where 


m i = f Q f 1 + “e 5 (0] yi 2 (l)<3-l 

p i( T ) = f 0 t 1 + “e y s (^T)de 

Sim ilar ly, substitution of expression ( 10 ) into the initial conditions 
(eq. (9)) and taking into account equations (12) lead to the following 
conditions : 


0i(O) 


P 3 (°) 

m i 


01 ( 0 ) 


Pj(0) 

“i 


(i = 0, 1, 2, . . .) 
(i = 0, 1, 2, . . .) 


( 20 ) 


A simple formula for the generalized-mass integral m^_ for 
i = 1, 2, . . . has been presented in reference 9 for uniform beams 
having any of the usual end conditions (free, pinned, or clamped) but 
without concentrated masses. In terms of the dimensionless quantities 
defined herein, the extension of this formula to beams with a concentrated 
mass 5 c at | = 0 is 

[l + “c 5(l)] y i 2 (l)<ii 

\ 

= £ 5b yi 2 (°) + - 2y 1 , (i) ^'"(i) + yi" 2 ( 1 )l 

* "* * 

(i = 1, 2, . . .) (21) 

For the present case where the end (| = 1 ) is free, equations ( 21 ) reduce 
to 

m ± = y ± 2 (0) + yi 2 (l)J (i = 1, 2, . . .) (22a) 

The rigid-body generalized mass (i = 0 ) is seen to be 
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“0 = C 1 + “c)^0 2 ^ 22b ^ 

Some reduction of the generalized-load integral P^(t) may also he 
accomplished in general terms for i = 1, 2, . . . . The quantity 

[l + m c 6(!)]y ± (0 may he replaced hy yj_^(6) (eq. (l5))j then, 

successive integrations hy parts and application of the boundary condi- 
tions (eqs. (l^) and (l6)) reduce the integral to 

p i( T ) = -^2 f yi(l) d l (i = 1, 2, . . .) 

kj. JO 

Substituting from equation (13) and recalling that, in natural vibration, 
the inertia loads on a free-free beam are self -equilibrating yield, 
finally, the generalized load 

P^Ct) - • f q(|j T ) yi(l) d 6 (i = 2 , • . .) ( 23 a) 

k i Jo 

On the other hand, for i = 0, the quantity Pq(t) is most simply- 
expressed as 

p 0 (t) = y 0 f ( 2 3b) 

Jo 

It might be noted here that, in the usual method of normal modes, the 
expressions for generalized force corresponding to equations (23a) do 
not have the factor l^k^ 2 . This Is one manifestation of the more rapid 

convergence of the Williams method. 

The problem now requires direct integration of equations (ll) and 
( 13 ) for the deflections y r and y g , solution of equations ( 19 ) for 

the generalized coordinates j)^, and solution of equation ( 15 ) for the 
natural modes of vibration y^, -with each function satisfying the desig- 
nated boundary or initial conditions. Direct integration of equation (ll) 
with the initial conditions (eqs. (12)) taken Into account yields 

y r ( T ) = 1 _ f f f q( !>T)dg(dT) 2 

1 + ^JqJqJo 


(2k) 
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Substituting equation (ll) into equation (13) and integrating four times, 
taking into account the boundary conditions on y g (eqs. (l4)), yields 

the following result: 


y s (i/r) 


r 6 f l r l r 5 ste.TKdD 1 

J 0 0 J 1 J 1 


1 (i k _ i 3 
1 + mc\24 "5” 


+ 



q.(l,T)d| 

( 25 ) 


The solution of equation (19) , satisfying also equations (20), is readily 
obtained by means of the Laplace transform. The result is 


0 ± (t) =-^- + ^/ 0 p i( 0 ) sln k i( T - 0 ) de (i = 0 , 1 , 2 , . . .) 


(2 6) 


Finally, the natural-mode shapes y^ and the corresponding frequency 

equation are derived in reference 10. These results, including the 
natural-mode shapes and the frequency equation, are summarized for easy 
reference in table l(a) . 

Relations are also given in table I (a) for the moment M(|,t) and 
shear V(|,t) obtained by substitution of the deflection response into 
equations ( 5 ) and ( 6 ). (It is also possible to obtain these quantities 
by Integrating the total load as 



However, some care must be exercised in using these formulas when the 
load function is discontinuous in time or has discontinuous first deriv- 
atives with respect to time . ) 

Other conf igurations . - The response of a uniform free-free beam 
without a concentrated mass is given by the results in table l(a) with 
m c = 0. The response of a cantilever beam may also be obtained from the 

response of the free-free beam with the concentrated mass by a limiting 
process in which the mass m c approaches infinity. Results for the 

cantilever beam are summarized in table l(b). For completeness, the 
Williams solution for a simply supported beam is shown in table l(c). 




# 


V 


t 
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Time -dependent boundary conditions .- It is worthwhile to point out 
that the method outlined in this report is directly applicable to the 
solution of problems with nonhomogene ous time -dependent boundary condi- 
tions . Such problems require the separation of the solution into two 
parts , one satisfying the time-dependent boundary conditions and the 
other capable of being expanded in terms of time -independent functions 
such as the natural modes of the beam. (See, for example, ref. 11.) In 
the Williams method, this separation is already made and time-dependent 
boundary displacements or forces are simply introduced into the boundary 
conditions imposed on y s or into the equations for rigid-body 
displacements . 

Consider, for example, a uniform beam fixed at one end and given a 
variable displacement at the other, such that its differential equation 
and boundary conditions are 

y^CtjT) + y(| ,r ) = q.(| ,t) 
y(o,T) = y' (0 ,t) = y" (1 ,t) = o 


y(l,r) = s(t) 

The solution would be assumed in the form of equation (10) but with 
y^ = 0 since there is no rigid-body translational freedom in this case. 
The static portion of the solution would be determined to satisfy 

_ IV, 

y B (|,t) = q(|,T) 

and 

y s (0,T) = y g '(0,T) = y s "(l,T ) = 0 
y s (l,T) = s(t) 

while the expansion functions y^ (where i = 1, 2, . . .) are the solu- 
tions of 

y^Cl) = 
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and 


y ± (o) = ^'(0) = y ± "(i) = 54(1) = o 

(the natural modes of a clamped-pinned beam) . In order to complete the 
solution, the generalized coordinates corresponding to a beam initially 
at rest and unstressed would have the usual form 

$i( T ) Pj_(9)sin k ± (r - 9)d6 (i ■ 1, 2, . . .) 

where 


Pi(t) = -iff r yi(|)d| 

J 0 

=4. = f Q = “^2 ^i^ 1 ) yi"^ 1 ) 

Similarly, a uniform free-free beam with a specified time-dependent 
displacement g(-r) at its center moves according to 

y^T) + = q(|,T) 


y(0,r) = s(t) 


y* (0,T) = y" (1 ,t) = y ,M (l,T) = o 


In this case 


y r (*r) - s(t) 


and y a (i,*r) is determined frcm 
s 


y 6 IV a,r) - q(l, t) - g(T) 


and 

y s (0,T) = y s '(0,T) = y s "(l,r) = y s '"(l,T) = 0 


V 


r 
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while the natural modes y^ (where 1 = 1, 2, . . .) are the modes of 

a cantilever beam. Or if, Instead of y(0,T) = g(T), there is given the 
force boundary condition y'" (0,r) = h(r), the rigid-body motion is deter- 
mined from 



The static solution y s is taken to satisfy 

IV 

y s (i > T ) = «[(S/r) - y r (T) 

y B (0,r) = y g ’ (0,t) = y B ”(l,T) = y s ,n (1,t) = 0 

and the modes y^ (where i = 0, 1, 2, . . . ) are the natural modes of 
a free-free beam. In- this case, it can be shown by integrating the dif- 
ferential equation governing y g that y s ,,, (0,T) = 1i(t) . 

Thus, the treatment of problems with time-dependent boundary condi- 
tions involves no special separate procedure when the Williams method Is 
used. 


Timoshenko’s Beam Theory 

Basic equations .- In the elementary beam theory, deflection occurs 
only by virtue of the rotation of the beam elements and only their trans- 
lational inertia is taken into account. The Timoshenko beam theory 
(ref. 9) permits additional deflection due to transverse shear and accounts 
also for the rotational Inertia of the beam elements. According to this 
theory, the motion of a beam subjected to an applied load of Intensity 
q(x,t) is governed by the equations (see, for example, ref. 6): 


— El ^ + A s g(— - t] 
dx bx. \Bx / 


2 n 
mr — 7T = 0 




dx 


A s o(k 
s Vcix 




(27) 


where t}t is the rotation of the cross section (see fig. l), r is the 
radius of gyration of the cross section, and A S G is the stiffness in 
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transverse shear. The effective shear-carrying area A s differs from 

the total area because the shear stress is not constant over the cross 
section. The bending moment M and transverse shear force V are 
given by 


M = -El (28) 

and 

v = AsG (i^ ■ *) (29) 

For uniform beams, these equations may be written in the dimension- 
less forms 



2 ~ 

^ t = 0 

> 

-q. 

- 


m = --r 


(50) 


(51) 


V = -^p(y' - t) 
k 2 


(52) 


where the transverse shear coefficient k a = is a measure of the 

3 1 \Ja 8 g 

freedom of the beam to deflect through transverse shearing action and the 
rotary inertia coefficient k^j = j is a measure of the rotational 
inertia per unit length. 


Note that the functions y and are both necessary for adequate 
definition of the deformation of the beam. Since these go hand-in-hand, 
the terms "solution" and "response," as used herein, will apply to these 
functions collectively and the single notation y(l.> T ) j 1 Kl>' r ) will be 
used to specify both functions . 


Symmetrical free-free beam with concentrated mass .- For the applica- 
tion of Timoshenko's theory to the symmetrical motion of a uniform free- 
free beam with a mass 2m c at the center | =0, attention will again be 
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restricted to the portion 0 ^ 4 1. As in the case of the elementary 

theory , the effect of the mass may he introduced into the differential 
equations if desired. However, for illustrative purposes, the differ- 
ential equations ( 30 ) will be left unchanged and the mass will be intro- 
duced in the boundary conditions; the boundary conditions then become 


*(0,t) = 0 
= 0 

y* (i> T ) - t(i,T) = 0 

-^2 y* (0,T) = m c y{0,r) 
k s 


(33) 


« 


\ 


Note that the location of the concentrated mass at 4=0 and the 
restriction to symmetrical motion exclude any effect of the rotational 
inertia of the concentrated mass. 


Hie beam is assumed initially at rest and undeflected; hence. 


y(S,o) = 0 

Hi, 0 ) .. 0 ^ 

♦( 6 , 0 ) = 0 

*( 4 , 0 ) = 0 ^ 


(34) 


With the problem thus completely defined by the differential equations 
(eqs. 30 ), the boundary conditions (eqs. (33)), and the initial conditions 
(eqs. (34)), the solution may be obtained as follows. 


Assume that 


y(fc, T ) 


*(l/0 


y r ( T ) + y s (s> T ) + Y ^(t) y^i) 

i=6 


♦ B (S/0 + ^T 0 ± (t) ti(l) 
i=0 


( 35 ) 
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where y r (r) is again the rigid-body translation of the beam, 

y s (l^) j^ s (£,t) is the static solution, and yi(|);ti(l) (where 

i = 0, 1, 2, . . .) are the natural vibration inodes. The rigid-body 
translation of the beam y r is governed again by the differential 

equation 


l - rh: f 1 

r 1 + mo Jo 


and the initial conditions 

y r (°) = y r (o) = 0 

Hie static solution is determined to satisfy 


V + 

s k, 


“ %) “ 0 
s' ' 

h$s' - *s)' = “5 + 


and the cantilever boundary conditions 


•"V 


y s (0,r) = 0 
* S (°> T ) = 0 
V(1,t) = 0 
y a * (1 ,t) - * b (1,t) = 0 

The mode shapes 3^(1 ) satisfy 

V + ri(h' - *i) * h % l \ 

K s 

Afo.' " *i) + k i 2 ^i = 0 


= 0 


(36) 


(37) 


(38) 


(39) 


(40) 


* 


i 


v 


* 
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the boundary conditions 


* ± (0) = 0 
V(i) = 0 
yi'Ci) - tid) = o 
^2 = ^ 1 ( 0 ) 


m 


and the orthogonality relation 


r* • 

f Q ' [ x + a c yj(s) + kRi^Ci) *j(t) 


d| = O 


(1 t J) 
(te) 


The derivation of this orthogonality relation is shown in the appendix 
along with the solutions to equations (40) and (4l) . 

Substituting equations (35) into the differential equations (eqs. (30)) 
and utilizing equations ( 38 ) and (40) reduce the differential equations to 


t (?i + H 2 P±)h - -y B 


i=0 


(*3a) 


CO 

£ (?i + VW)*! = -t E 


(^5b) 


i=0 


Multiplying equation (b3a) hy jjL + m c 8 ( £ )J y^ and equation (43b) by 

, adding the two equations, and integrating the sum over the range 
0 £ < 1 yield the result 


h + *±\ 


(i = 0, 1, 2, . . .) (kk) 


v 
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which takes into account the orthogonality relation (eq. (42)). The 
generalized mass and generalized load appearing in equation (44) are 


m. = 




4 


and 


1 

p i( T ) - f t 1 + “c B(£)]yiy 8 + kRi^i^g 

d 0 


<4 


respectively. 

By a similar process the initial conditions (eqs. (34)) become 


i(0) —^21 (i = 0, 1, 2, . . .) 


0 ± (O) = - 


m i 

Fj(0) 

m i 


(i - 0, 1, 2, . . .) 


(45) 


It is shown in the appendix that the generalized masses mj_ of the 
given free-free beam can be evaluated for i = 1, 2, . . . as follows: 


“i 




(i “ 1> 2, . . .) (46a) 


The remaining generalized mass mQ reduces, as in the elementary theory, 
to 


2 

n>Q = (1 + 


(46b) 


since the symmetric rigid-body mode has no rotational component Vq. 
Further, the generalized force P^(t) for i = 1, 2, . . . may be 

reduced, by a process of substitution from equations (38) and (4o) and 
integration by parts, to 


?i( T ) = -^2 yiCl) i(l> T )<4 


(i = 1, 2, . . .) (47a) 
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and the rigid-body generalized force P Q is 


* 0 ( T ) 




(4Tb) 


On the basis of the assumed form of the solution expressed in equa- 
tions (35), the problem of determining y(S ,t) ;*(£ ,t) has been replaced 
by a number of component problems requiring determination of the func- 
tions y r (r), y s (i,T)^ s (|,r), 0 ± (t), and y^OitjU). The solutions 

of these component problems must now be obtained. 

For y r (r), integration of equation ( 36 ) in conjunction with equa- 
tions ( 37 ) yields 


y (t) = — f f f *(l> T )d|(dr) 

r i + DjJoJoJo 


(48) 


For y (|,T);t s (^ T )^ substituting y r fr cm equation ( 36 ) into equa- 
tions ( 38 ) and integrating, in conjunction with the boundary conditions 
(eqs . ( 39 ) ), gives 


a ¥-s)-fe- T + t)] / 0 X 5(, ’ T,dt 


l + m c 




(49) 


Since equations (44) and (45) are identical to equations (19) and (20) of 
the elementary solution, the generalized coordinates <f> ± (r) are again 

given by equation ( 26 ) . Ihe solution of equations (40) for the natural 
modes y-j_(l ) ;%(£ ) is given in the appendix. 

Hie solution just obtained and corresponding solutions for canti- 
lever and simply supported beams are summarized in table XX. 
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TWO NUMERICAL METHODS OP SOLUTION 


Traveling-Wave Method 

A traveling -wave method for calculating the response of a structure 
to transient loads can he devised only if the motion of the structure is 
governed by differential equations of the hyperbolic type. The simplest 
beam theory which completely fulfills this requirement is Timoshenko's 
theory, which includes the effects of both transverse shear and rotary 
inertia . 


In developing a traveling-wave method of solution, it is convenient 
to first replace the Timoshenko partial differential equations (eqs. (27)) 
with the following four equivalent ordinary differential equations written 


along 

Along 

four 

1+ 

characteristic lines I+, I-, H+, and H- in the 

where ^ dM + mr 2 dil - V dt = 0 

dx c-j_ c^ 

x,t plane: 
(50a) 

Along 

I- 

where 

dt _ l . 
dx c-,/ 

— dM - mr 2 dfi + V dt = 0 
C 1 

( 50 b) 

Along 

11+ 

where 

dt = _1_. 
dx c 2' 

dV - m dv + ^mCgft + q)dt = 0 

( 50 c) 

Along 

II- 

where 

dt __1_. 
dx c 2* 

^ dV + m dv + ^mCgO - q^dt = 0 

(50d) 


The derivation of equations (50) is given in reference 6. The dependent 
variables are the moment M, the shear V, and the linear and angular 


dy dijf 

velocities v = — and f2 = — - . 

St St 


The quantities 


C 1 = 



and 


= JM 

y m 


are the propagation velocities of discontinuities in moment and shear, 
respectively (phase velocities of disturbances with infinitesimal wave 
length) . In each equation, the total differentials specify infinitesimal 
differences in the designated characteristic direction. 


For any given beam, the slopes of the characteristic lines are known; 
hence, closely spaced networks of characteristic lines may be drawn in the 
space-time plane. Various schemes for the approximate step-by-step inte- 
gration of equations (50) over such networks are possible. In general, 
all require some form of interpolation since Timoshenko's equations have 
two characteristic nets. (The particular case where the two nets coincide, 
C 1 = c 2> ^as keen treated in detail for uniform beams in ref . 6 . ) One 
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integration scheme is briefly described in this section. Attention is 
restricted to a uniform beam for which the characteristics are straight 
lines . 

The beam is divided arbitrarily into segments Ax (see fig. 2 ), and 
the time interval is taken according to At = -i- Ax. This specifies a 

lattice of points in the space-time plane at the intersection of the 1+ 
and I- characteristic lines. Consider a general point 1 (fig. 2) from 
which characteristics of both families have been drawn backwards in time. 
The 11+ and II- lines have steeper slopes ^ since c 2 < cj) and termi- 
nate at points 2' and 4' . Then the differentials in equations (50) may be 
replaced by the appropriate finite differences and the following equations 
are obtained: 


- Mg) + mr 2 ^ - n 2 ) - ^( Vl + V 2 ) = 0 (51a) 

i(Mi - Mi,) - mr 2 (0 1 - 0^) + + V^) = 0 (51b) 


- v) 

| - m( Vl - v 2 .) 

+ mc 2 ir( ft l + ft 2’> 

1 + (oLp + 

0 

u 

(5lc) 

4( v i - M 

+ m^ - Vi,,) 

+ mc 2 ££(% + n v ) 

- f (*1 - 

qi,») = 0 

(5ia) 


It is assumed that M, V, v, and & are known at points 2 , 3> 
and it-. Parabolic interpolation formulas are substituted Into equa- 
tions (51c) and (51d) to give the quantities at points 2 ' and V in terms 
of their values at points 2 , 1 5 , and 4. Then equations (51) become four 
equations for the four unknown quantities V;]_, v^_, and In 

terms of known values of M, V, v, and fl at each of the lattice 
points 2 , 3 } and 4. These may be solved to obtain a matrix recurrence 
formula; however, a simplification may be introduced based on the fact 
that quantities at points 2 and 4 have already been determined to satisfy 
the following equations: 


i(M 2 - Mj) - »r £ (n 2 - n,) + ^(v 2 + v 5 ) - 0 
bx( M 4 - *b) + - <b) - ir( v 4 + v 5 ) = 0 ^ 


(52) 
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Combining equations (52) with equations (51a) and (51b) leads to the 
result 


mr 


! (°1 -%-% + %)- f (V! - v 5 ) = 0 


(53) 


where M has been eliminated. Now equations (51c), (51d), and (55) con- 
stitute three equations for the three unknowns ft^, V-j_, and in 

terms of known values of ft, V, and v at each of the points 2, 3, 4, 
and 5 . Solution of these equations yields the recurrence formula 


ft n 




&2 


ft ^ 


% 




h 

v 2 

v 2 

+ [ a 3] 


+ [ a 4 ] 

^4 

^4 

+ [ a 5 ] 

h 

^5 

+ [qJ 

*2 

*3 

*4 


(54) 


where the quantities have been made dimensionless = — = ni 2 |/,g- 




and v = = vl \ & \ and where 


St 


El 


St 



C 1 

_£2 K 
C 1 ^RI 

O 

IV) 

ro 

o 

ro 

ro 

H 

c2 

C 1 

^ ^RI 


If + K?) 




0 

c 2 

1 0 
0 1 + K 2 


* 


£ 
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2 - 





-if V* 1 + k 2 ) 


Co 

— K 

C 1 



KrjCI + K2) 


f2 JC_ 
C 1 k RI 



^(1 + K2) 

C 1 


M- 


-1 

c l 

0 


_£l 


K 


'2 

K 2 


0 

0 

0 


M 


“e/ & ^(l 4- K 2 ) 


“1^ 


— K 


" I *51 C X 
| ^I 2 % K (i + ^ 



K = 


1 c 2 At 


Note that the more appropriate parameter 

instead of k . 

s 




has been used here 


The response of a beam may now be obtained by the repeated applica- 
tion of equation (54) except that, as is indicated in figure 2, special 
formulas which take into account the particular boundary and initial con- 
ditions of the problem, must be derived for boundary and initial points. 
In addition, it must be remembered that the characteristic lines are pos- 
sible loci of discontinuities in the dependent variables or their deriva- 
tives. (See, for example, ref. 6.) Such discontinuities will arise, for 
example, if there are concentrated loads (or imposed velocities) that 
have histories which are discontinuous or have discontinuous derivatives. 
Discontinuities in M or D are propagated with velocity c^ (the 
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locus of such, a discontinuity is shown by the dashed characteristic lines 
in fig. 2 )} discontinuities in V and v propagate with velocity Cg. 

The magnitudes of discontinuities may be predetermined through the appli- 
cation of equations (50) in a manner which is illustrated in reference 6. 
Thus, in general, discontinuities may be added as they are encountered in 
the step-by-step' solution. In the scheme which led to equation (54), it 
will be noted (fig. 2) that special consideration is necessary for pointB 
just above the locus of a discontinuity. Special formulas are certainly 
required where there is a jump in one of the functions and it may be 
desirable to account also for discontinuous first derivatives. 

Once ft, V, and v have been determined at a point, M may be 
obtained at that point by integrating equation (50a) or equation (50b) 
along the proper characteristic from some boundary where M is known. 

It should be. mentioned that the selection of the 1+ and I- lines 
as the basic network is based on (rather intuitive) considerations of the 
stability of the numerical procedure with regard to propagation of errors. 
(See ref. 12.) It is assumed that the domain of dependence imposed by a 
numerical procedure should at least encompass the total theoretical 
domain of dependence. This would not be the case if the steeper 11+ 
and II- characteristic lines were utilized as a basic grid. 


Houbolt's Method 

A cantilever beam acted upon by a series of concentrated loads p 1 , 
P£, . . . p n . Is shown in the following sketch: 



0 


l y \[ 

12 5 

Station 



Such a beam has the deflection 


y c | = [ G ]M 


(55) 
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where the subscript c is used here to indicate that the deflections 
are cantilever deflections measured with respect to station 0 and where 
[G] is a matrix of stiffness influence coefficients. The inverse equation 
is 


IP I - M'Vcl (56) 

The deflections of a free-free beam may be expressed in terms of the 
cantilever influence coefficients [G]. For the free-free beam the deflec- 
tion is y = y 0 + y c where y 0 is the deflection at station zero (the 

center of the free-free beam); hence, equation ( 56 ) becomes 


|p| = M -1 |y I - yoM " 1 


(57) 


But, on the free-free beam, there is the additional load p Q . From the 
condition of overall equilibrium, p^ is given by 


n 


1 

1 


?0 --I Pi ■ " L.i 1 • • • iJO^J |y| + ybL 1 1 * • • -*J 03 


-1 


(58) 


i=l 


and equations ( 57 ) and ( 58 ) may be combined into a single matrix equa- 
tion for the loads at stations 0 to n of a free-free beam. This 
equation may be written 


|P| = Miy | 


(59) 


a (bj 

I b | Cg] _1 


where 
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1 

1 

M = -[ g ]" 1 • 

1 

Lbj = -Li i . . . ij [g ]" 1 

i 

i 

a = [l 1 • • • iJCG]" 1 • 

1 

ana where the vectors ] p | ana f y | now contain p Q ana y Q terms, 
respectively. 

Consiaer the motion of a free-free beam subjected to the applied 
distributed load q.(x,t). The equivalent concentrated loads at the 
stations 0, 1, ... n are, at any time, 

| P I = I>](kl - MJy|) ( 6 °) 

where £m] ls the diagonal mass matrix, and [w] is the weighting matrix 
(ref. 15): 
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7 6-1 
2 20 2 
2 20 2 



Ax 

2 k 


2 20 2 


2 20 2 
-I67 


The essential feature of the Eoubolt procedure (ref. 8) is the method of 
expressing the second time derivative jf. The acceleration at time 
t. = j At is written as follows: 


W i = ' 5ly|j ' 1 + 4|ylj - 2 " ly ' J 4 


(61) 


and is obtained by passing a third-degree curve through the points at 
t = tj, ^j-l* ^j- 2 > eac ^ 1 station x^. Substituting equa- 

tions (60) and (6l) into equation (59) (written for time t = tj) and 
solving for yj leads to the recurrence equation 


|y|j - Wklj + ^i^( 5|y| j-i ‘ 1,|yl J-2 + Wj- 5 ) 


(62) 


where 


-.-1 


[B] = 


[A] + 


(At)‘ 




[W] 


[c] = KJM 
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Thus, the Houbolt method, simultaneously determines all the deflections 
y(x^,tj) (where i = 0, 1, . . . n) in terms of the deflections 

y(x±,tj_i) , y(x ± ,t J-2 ), and y^x i ,tj_^^. A concentrated load 2f(t) 

at the center of the free-free beam may be included in equation ( 62 ) by 
adding to the right-hand side the term 

fjlal (63) 

where 


Id | = 



1 

0 

0 


0 


The solution of a problem by repeated application of equation (62) 
requires that the initial conditions of the problem at time t = 0 be 
expressed in terms of fictitious ordinates | yj t_^_ 

and t_ 2 - This is accomplished by expressing the time derivatives | y| Q 

and IyIq in terms of a third -degree curve passed through points at 

t = t^, 0, t ^ and t _2 at every station. For the case where the 

beam is initially at rest and the applied loads Increase continuously 
from zero initial values, the initial conditions 


i y ' 0 ‘ |Jl o = H ■ 0 

are converted by this process to 

|y| = 0 

0 

M_i “ -Mi ’ 

|y |_ 2 = - 8 lyl 1 


(64) 


(65) 
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Then, equation ( 62 ), applied for j = 1, yields 

|y| 1 = [B]|q| i-^ Cc]|yl i 

which may he solved for |y|^ to obtain 

-1 

|y|i = [El + H |aJ i (66) 

Bie application of equation ( 62 ) for j = 2, 3, . . . is now straight- 
forward. All that remains is the determination of the cantilever influ- 
ence coefficients £g]| . 

In this connection it should he pointed out that, although no 
restrictions have been made on the beam theory to he used, the applica- 
tion of the method, as formulated, with Timoshenko's theory requires 
that the deflections y he interpreted in a general sense and include 
also the rotations of the cross sections i|r. Thus, with two quantities 
to he determined at each station, the order of the matrices is 2n and 
the required computational labor is roughly four times that required with 
the use of the elementary beam theory. A compromise which affords 
increased accuracy over the elementary theory, yet avoids thiB .large 
increase in computational labor, is the use of a theory which contains 
transverse shear freedom hut no rotary inertia. With no associated 
inertia loading, the rotations i|r need not he explicitly included in 
the step-by-step dynamic analysis and do not appear in the recurrence 
formula, equation (62) . The cantilever influence coefficients are deter- 
mined as follows on the basis of this latter theory. 

The influence function (Green's function) G^xjx^J is the solution y 
of the equations 



r (6 t) 



and the boundary conditions 
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y(0) = ♦(£>) = |*(l) = ^(2) - +(2) = o 


(68) 


However, the deflection may he written as the sum of bending and shear 

dy-n 

contributions, y = y B + y g with i|r = ^-=-,-and it is expedient to write 
for the influence function 

G(xjxi) = Ggfxjx-j.) + Ggfxjx-J ( 69 ) 

where Gg and Gg are the solutions of the differential equations 
(equivalent to eqs. (67)) 


h A z G h G s - -s(x-xi) 


(70) 


8x s 8x 

and boundary conditions (equivalent to eqs. ( 68 )) 

°bM ■ - Jr M ■ h EI ■ 0 


°s(°* x i) ■ ^(^i) ■ 0 


(71) 


For given distributions of bending and shear stiffnesses, equations ( 70 ) 
may be integrated directly in conjunction with the boundary conditions 
(eqs. (71)). For a uniform beam, the resulting total influence functions 
are: 


rfxjx-, ) = x + -^r 
\ ’ AG EI 

s 


EI \ 2 






«s _i_ x + I L 
A S G 1 EI\ 2 


*x 




(x < xj 
( x > x l) 


(72) 
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In dimensionless terms, equations (72) may be written 

3 (‘> s i) " k s 2 s + 4(h - 1) (* < «i) 

* k s 2| i + %■(? ■ t) ( 5 > 6 a) 

and, for a uniform beam, the recurrence formula, equation ( 62 ), 




a |j)J 

1*1 B8 




‘ (73) 

* 

becomes 

(7*0 
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7 6-1 
2 20 2 
2 20 2 



2 20 2 
2 20 2 
-16 7 


The elements G-j_j of the matrix of influence coefficients £gJ are seen 
to he 

hi - k s 2i *s + y( j • 3) (as)3 (i < J) 

- k g 2 J A5 + £(1 - |)( AS) 5 (1 > j) 

where i designates the row and j , the column. 

RESULTS AND DISCUSSION 


Modal Solutions 

For illustrative purposes, example computations have teen made for 
the case of a uniform free -free team, for which =0.1, k g =0.2, 

and m„ = 0. The team is subjected to the applied concentrated load 

q(|,T) = 2?(t)5(|) 

where f(T) has each of the time variations shown in figure 3* The 
given values of k^j. and k g are appropriate to a team having a solid 
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rectangular cross section and. a ratio of half-length l to depth, equal 
to 2 . 887 . Hie calculated response has "been limited to the history of 

transverse shear at the point g = Moment calculations are omitted 

because they do not provide as severe a test of the analytical methods. 

The point | = ^ is chosen arbitrarily, since the location corresponding 

to maximum transverse shear is not known in advance. 

The responses to the step and ramp-platform loads were obtained on 
the basis of elementary theory from table l(a) and on the basis of 
Timoshenko's theory from table Il(a) . In each case, six modes were used 
in the expansions. The responses to the other three functions (figs. 3(c) , 
3(d), and 3(e)) were obtained by superposition of the step and ramp- 
platform results. The resulting shear histories are shown in figures 
5 , 6 , 7 , and 8 up to a time corresponding approximately to the period of 
the first natural mode of vibration of the beam. 

The long dashed curve in each figure is the static portion of the 
response. Thus, the largest dynamic overshoot factor, achieved by the 
step loading (fig. 4), is approximately 2 . 7 - 

It will be noted that the blast pulse load (fig. 3(e)) has only one- 
half as much impulse as that contained in the square and triangular pulse 
loads (figs. 3(c) and 3 (<3L) ) * Hence, if the responses to the three pulse 
loads (see figs. 6 , 7 , and 8 ) are to be compared on the basis of equal 
input impulse, the response to the blast pulse must be doubled. On this 

basis, the blast pulse is seen to cause the highest shear stress at | = 

An indication of the convergence of the modal results in figures k- 
to 8 is provided by the bar graphs in figure 9- On each graph, the heights 
of the bars correspond to the magnitude of the static portion of the 
response (zero frequency) and to the amplitudes of the terms in the series 
expansion for the dynamic portion. (For each load, r is sufficiently 
large so that the load function has attained its constant value . ) The 
bars, thus, represent the maximum possible contribution of each term to 
the total. Note that, for the loads of long duration (the step and ramp- 
platform functions ) , the static part and the first term contribute a pro- 
portionately large share of the response and sufficient accuracy could be 
obtained with only three modes. On the other hand, the adverse effect of 
reducing the load duration is illustrated by the bar graphs for the 
responses to the three pulse loads. In each case, no convergence is appar- 
ent for the first few modes and, in the cases of the square and blast 
pulses, there is some doubt as to the adequacy of the even six modes, 
particularly with the use of the elementary theory. 


r— I |CV1 
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Further evidence of the effect of load duration on convergence and, 
in addition, an indication of the effect of load distribution are given 
by the following cases of the response of a uniform simply supported beam 
to various loads: 

Case (l): For a uniformly distributed step load q.(|,T) = _^(t), 

00 

M(&,t) =-l£ + l- iL Y i. sin iit| cos i 2 :! 2 ^ 2 

2 2 sP L ±5 

i=l,3i • • • 

oo 

V(£,t) = -t + \ " ^2 ^ JS cos ijt 6 cos i 2 ’! 2 !' 2 

i=l,3,... 


Case (2) : For a step load concentrated at the center 

= 5 (t-0/( T ), 

CO 

M(6,t) = | - ^ ^ (- 1 ) ^ sln i3t 6 cos i 2 ^ 

i=l,3, • • • 


»(t,T) £ (-D 


1-1 

2 


i cos iit£ cos i 2 it 2 T 


i=l,3. 


Case (3): For a uniformly distributed impulse load <3.(£,t) = 6(t), 


S1(|, t ) = i Y j sin iit| sin i 2 jt 2 T 

i = i , Y , . . • 


V(6,r) = k 


oo 

I 


2 2 

cos irt§ sin i it t 


i=l,3. 
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Case (4) : For an impulse load concentrated at the center 

q.(£,T) = ”i) 8 (r)' 

00 i_l 

M(g,r) =2 (-1) 2 sin lit I sin i 2 rt 2 r 


1 = 1 , 3 , 


V(6,t) = 2rt £ ( 

i=l,3, • • • 


i-1 

-1) i cos irt| sin i 2 rt 2 T 


These results are based on elementary beam theory. (The last two cases 
are simply time derivatives of the first two cases . ) The effect of load 
duration is illustrated by comparing cases (l) and (3) (or cases (2) 
and (4) ) . Changing from a load of infinite duration (the step load) to 
a load of zero duration (the impulse) introduces a factor i^ and, hence, 
reduces the rate of convergence (and, in fact, produces divergent series 
in both case (3) and case (4)). Similarly, the effect of Bpatial dis- 
tribution may be seen by comparing cases (l) and (2) (or cases (3) and 
(4)). Changing from a distributed to a concentrated load introduces a 
factor i and hence reduces the rate of convergence. The apparent change 
in the sign of half the terms is not significant since each series is 
essentially an irregularly alternating series (except at certain specific 
combinations of | and t) . 


The Need for Timoshenko ' s Theory 

Thg question of which theory should be used to determine the response 
of a beam to a transient load is intimately related to the convergence of 
the result. This is because the secondary effects of transverse shear and 
rotary inertia become increasingly important for the higher modes. Note 
in figure 9, for example, the growing disparity between the natural fre- 
quencies of a uniform free-free beam calculated on the basis of the ele- 
mentary and Timoshenko theories. Thus, if it is determined that, for a 
given beam subjected to a certain load, modes strongly affected by trans- 
verse shear and rotary inertia contribute a large share of the response , 
it is unlikely that elementary theory will yield correct results. In case 
the given beam is a complicated nonuniform structure, a rational procedure 
for dete rmin ing the proper theory would be to consider a uniform approxi- 
mation to the given beam, quickly obtain the response of the uniform beam 
to the given load by elementary theory (table I) and investigate the con- 
vergence of the response, and at the same time to consider the influence 
of rotary inertia and transverse shear on the modes (as manifested by the 
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differences in natural frequencies obtained with the elementary and 
Timoshenko theories). 

This reasoning is generally, though not conclusively, confirmed by 
the results in figures 4 to 9* Only the first mode of the free-free beam 
shows good agreement between the frequencies as given by the elementary 
and Timoshenko theories. (See fig. 9*) Thus, the responses to the step 
and ramp-platform loads (figs. 4 and 5) , which depend heavily on the 
first mode and the static contribution, also show good agreement between 
the two theories. The responses to the pulse loads (figs. 6, 7> and 8), 
obtained with the two different theories, bear little or no resemblance 
to each other since they depend heavily on the higher modes. However, 
except for the response to the square pulse load (fig. 6), the two 
theories do yield about the same peak stress, (it is felt that the posi- 
tive peak achieved in the first half period should be given more weight 
than the negative peak achieved later, since the latter would be con- 
siderably diminished by material damping which has not been included in 
this analysis.) 

Boundary conditions also influence the need for a more refined 
theory. For example, although the elementary theory is adequate for 
obtaining the shear due to a step load on the free-free beam of figure 4, 
it cannot be used to obtain the shear at the center of the same beam 
where the input is a prescribed "step-velocity" of the point £ = 0 (the 
so-called problem of the "instantaneous arrest of the root of a moving 
cantilever beam"). In the latter case, elementary theory yields a diver- 
gent series (ref. l4), whereas Timoshenko's theory yields a modal solution 
which converges to finite values (ref. 6). 


Numerical Solutions 

The two numerical procedures, the traveling -wave method and Houbolt's 
method, have also been used to calculate the transient response of the 
uniform free-free beam considered in the previous sections to an applied 
ramp -platform load. In the calculations by both procedures, the beam was 
divided Into six segments (A| = O.I667). Results are shown in figure 5* 

For the traveling-wave method, the time interval is 


At = kjy- A£ = 0.01667 

Relatively simple boundary formulas (based on linear interpolation) were 
used in this case and discontinuities in the derivatives of V and v, 
arising from the discontinuities in the slope of the ramp-platform func- 
tion, were ignored. 
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In contrast with the traveling-wave method, the Houbolt method 
imposes no inh erent restriction on the selection of the time interval At 
in relation to the space interval A£. The time interval may he taken as 
large as is consistent with the desired accuracy. This freedom has been 
utilized in that the calcinations by the Houbolt method have been made 
with the time interval At = 0.03333* which is twice the time interval 
used with the traveling-wave method. In addition, the computations by 
the Houbolt method do not include the effects of rotary inertia and were 
stopped at a point just beyond the peak load. 

It will be noted that the numerical results in figure 5 approximate 
the (essentially converged) Timoshenko modal solution fairly well. Both 
numerical methods underestimate the peak stress, the traveling- wave method 
by 14- percent and the Houbolt method by 9 percent. 

The greater accuracy of the Houbolt result Is particularly significant 
since it was obtained with less computational labor due to the use of the 
larger time interval. It must be kept in mind, however, that economical 
use of the Houbolt method requires that it be applied in connection with 
the elementary theory or with the addition of transverse shear alone. If 
both rotary Inertia and transverse shear must be included, the labor 
required in applying the Houbolt method is quadrupled. Fortunately, 
rotary inertia is negligible in many problems. (See, for example, ref. 15.) 

It should be mentioned that the traveling-wave procedure, as so far 
conceived, has proved somewhat sensitive to minor changes in the scheme 
used to obtain recurrence formulas. (Only one scheme, the most success- 
ful so far found, has been described in this report.) A possible reason 
for this stems from the fact that the slope of the basic network of char- 
acteristic lines is dependent on rotary inertia. In fact, rotary inertia 
is necessary to give the beam equations the wave character essential in 
the conception of a traveling-wave method. Thus, in view of the relative 
negligibility of rotary inertia for many practical problems, this sensi- 
tivity is perhaps not surprising. In general, it must be concluded that 
a traveling-wave numerical method of analysis which is superior to the 
Houbolt method has not yet been devised. 

Numerical results have also been obtained for the response to the 
triangular pulse load by superposition of the ramp-platform results. 

These results are shown in figure 7* The results indicate that more 
degrees of freedom must be taken with both procedures to predict ade- 
quately the response of the beam to the given triangular pulse load, and 
from this example It appears that the modal method of solution Is to be 
preferred. However, the simple problems discussed herein do not portray 
the main advantages of numerical methods. For example, numerical methods 
are readily extended to apply to nonuniform beams and conveniently adapted 
to the use of modern high-speed computers . A fundamental characteristic 
of numerical methods is the replacement of professional engineering time 



ko 


NACA TN k2kk 


by routine computing time. Hence, numerical procedures are not to be 
condemned on the basis of the results in figure 7* The selection of 
the best method requires the consideration of all these factors in rela- 
tion to the specific problem. 


CONCLUDING REMARKS 


Williams type modal solutions, based on both the elementary and 
Timoshenko beam theories, have been given for the response of several 
uniform beams to a general transient load. The response to any specific 
load may be obtained from these solutions by performing a series of indi- 
cated direct integrations of the load function. Typical computed results 
have been shown for the shear response of a free-free beam to various 
concentrated loads. 

The convergence of modal solutions is shown to depend both on the 
history and distribution of the load. Decreasing either the duration 
of the loading or the region over which the load is applied reduces the 
rate of convergence and may produce divergence. 

The need for a more refined theory, as compared to elementary theory, 
is intimately related to the rate of convergence of the modal solution. 

If modes which are strongly dependent on transverse shear and rotary 
inertia contribute a large portion of the response, Timoshenko's theory 
must be used. 

Comparison of the Houbolt and traveling-wave numerical methods indi- 
cates that the Houbolt, procedure has many advantages over the traveling- 
wave procedure as so far conceived. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., February 5, 1958. 
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APPENDIX 

SYMMETRICAL NATURAL VIBRATION QF A UNIFORM FREE-FREE 
TIMOSHENKO BEAM WITH A CONCENTRATED MASS 


Natural Modes and Frequencies 

The differential equations and boundary conditions governing sym- 
metrical natural vibration of a free-free beam with a concentrated mass 
at its center may be written in the following dimensionless forms: 


t" + ~g(y' - t) + = 0 

-^(y' - *) ' + k 2 ^ = o 

k s 


(Al) 


*( 0 ) = 0 


(A2a) 


t'(l) = 0 


(A2b) 


y* (i) - *(i) = o 

y' (0) + m c k 2 y(0) = 0 


(A2c) 

(A2d) 


where y(|);^(5) is the natural mode and k is proportional to the 
circular frequency of vibration. Each of the solutions of equations (Al) 
has the form 


y(l) = Ae^ 

f(£) = Be 7 ^ 


where A and B are arbitrary constants. Substituting this form into 
equations (Al) leads to the biquadratic equation 
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* 4 + + '‘ex 2 ) 7 ' 2 + * 2 ( A e\l 2 - l) - 


(A3) 


and to the following relationship between A and B:. 

. i,2v 2 


B=£±ai. 

A 


Equation (A3) has the four solutions X = ±a and X = tip where 

■\ 


(A4) 


a = kl/| 


-{ k 2 + ^i 2 ) + - *ta 2 ) 2 - 4 


e = k # 


( k / + ^ I 2 ) + |/( k S 2 - k EX 2 ) 2 + 4 


(A5) 


The general solution of equations (Al) may be written in the form 

**< 

y(l) = c i cosh a 4 + Cg sinh a| + Cj cos pi- + sin p| 

a 2 + 2 / » 
iKt) = — (Cg cosh a| + C-^ sinh af-J + 

P 2 - k 2 k 2 ( \ 

§-(0^ cos P| - C 5 sin p|J 


(A 6 ) 


where C-^ Cg, C-j, and are arbitrary functions of k which must 

be determined in order to satisfy the boundary conditions (eqs. (A2) ) • 

Substituting equations (A 6 ) into equations (A2) yields four homo- 
geneous algebraic equations for C 2 , Cj, and C^. The existence 

of nontrivial solutions of these equations requires the vanishing of the 
determinant of coefficients . This criterion yields the frequency equation 


F(k) = | 7 sin p cosh a + cos p sinh a + 
P 


H In 


g - — g S - ^2 + ^2 _ £^Bin P sinh a + ^7 + i^cos p cosh aj = 0 


(AT) 
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where 


a 2 + k 2 k g 2 
7 p 2 - k 2 k s 2 

Hie solutions of equation (A7) are the natural frequencies of vibra- 
tion k^ where i = 0, 1, 2, . . . . For k = k^ the homogeneous 

algebraic equations may be solved for the relative magnitudes of the 
quantities C^, Cq, C^, and Cij.. For i = 1, 2, . . . , the resulting 

vibration mode shapes may be written in the form 


y^i) 


sin 


sinh ctj_ 


■ =• cosh a. £ - — cos p., | - 

Pj. 1 a i 1 


P, 2 -i,V, \T 

1 — ^cosh + 7 1 cos p^J cos p^ - cosh a^g + 

+ P-L L 


K ± (sin p ± | - ^ sinh a ± g^ 


(i = 1, 2, . . .) 


t ± (S) 


„ c £i! k sin Pj_ 

% 1®i 1 Pi 


P-I sin Pj sinh a* 

— 7 . sinh Oi g + — sin p., g - 

°i 1 Pi ^ “1 1 


_ P ^ 2 - k ^ 2 !- 2 / \ f / \ 

~ — = — — Icosh + 7 ^ cos p^J A^Icos p^g - cosh cc^gl - 

oj 2 + p-j 2 ' ' L ' 


^sin p ± 


Z + ^ 7 ± S±nh a i^ 


(i - 1 , 2 , . . 0 


where 


Ai = 


(3-? 

sin Pi - sinh 
cos p^ + — cosh a± 


Hie rigid body mode, corresponding to kQ = 0, has the components y^ (|) 
and t 0 C|) = 0 . 
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Orthogonality of the Natural Modes 

The differential equations (Al) are satisfied by any of the infinite 
number of natural modes and corresponding frequencies. Thus, for the 
ith mode. 


V + - 1) + - o 




♦i)’ + *1%. " 0 


Let the first of these equations be multiplied by the rotational com- 
ponent of the jth mode and the second by the translational com- 

ponent y^ . If the resulting equations are added and integrated over 
the beam length, there results 



Integrating by parts the first two integrals on the right-hand side of 
the preceding equation yields: 




♦i)*j + V+j 


/q 1 V V ( A 9) 


This process is valid also if the roles of the ith and jth modes are 
reversed. Interchanging i and j in equation (A9) and subtracting 
the result from equation (A9) leads to 
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( k i 2 - k j 2 )/ 0 1 (vj +k BiVj >e = h)h 


(A10) 


If now the boundary conditions (eqs. (A2)) are imposed, equation (A10) 
Is found to reduce to 

( k ± 2 - k a 2 ) J 0 (hh + = -( k i 2 - k /H ? i (0) ? j (0) 


or 


f Q j[j- + £ c 5 (£)Jyi(5) yj(U + 4 ri 2 %(i) ^(6) 


d| = 0 


U i i) 

(All) 


since 


5c 5(|) y ± (g) y^g) d| = m c y^O) yj(O) 


Equation (All) is the orthogonality condition satisfied by the natural 
vibration modes of a uniform free-free Timoshenko beam with a concentrated 
mass at the center. 


Determination of the Generalized Mass 
The determination of the generalized mass 


^ ■ r 


[l + i c 5(6)] y^U) + kg/ */(!> 


dl 


by direct integration is a somewhat laborious process for i -f 0. Fortu- 
nately, m^_ can be expressed in terms of certain boundary values of the 
mode shapes by the application of a limiting process to equation (A10) in 
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which the functions y^ and. are considered as continuous functions 

of k. Eius, if 

kj = kj_ + dk 

and 

?j = ?! + ay = f ± + g dk 
“ +i + d* = tj. + dk 

equation (A10) becomes, in the limii; as dk approaches zero. 



(A12) 


This equation is applicable to uniform beams and may be extended, if 
desired, to nonuniform beams. On substitution of the boundary conditions 
(eqs. (A2)), equation (A12) reduces to 

f 0 [4 2 + VV] d 8 - - s c 4 2 <°> - 


1 

2k ± 



(i — 1, 2, . • .) 


Hence, the generalized mass is given by 




(i - 1, 2, . . .) 


(A13) 
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Note that only the second boundary condition (eq. (A2b)) is altered 
by differentiation with respect to k. This arises from the fact that 
only this boundary condition depends for its satisfaction on the frequency 
equation (A7)j that is, the mode shapes (eqs. (A8) ) satisfy the other 
boundary conditions for any value of k but satisfy equation (A2b) only 
for k = sirfce 

t'(l) = C {f - - k 2 k s 2 )F(k) 


where F(k) is defined by equation (A7). It is only by virtue of this 
dependence of one or more boundary conditions on the frequency equation 
that equations (A12) yield a value of mi. Hence, it must be concluded 

that, for a beam for which none of the boundary conditions depend for 
their satisfaction on the frequency equation (as, for example, a simply 
supported beam), equations (A12) are not applicable. For such a beam, 
mi is determined by direct integration. 
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TABLE X— RESPONSE at A lattPCBM EIEHEKEART BEAM TO A CEHERAL LOAD 


(a) Symmetrical free-free 'bean vith a concentrated Base 













Analytical expression 


y r (T) + y s (t,T) + J ^ t (r) j x (|> 

i-0 


m 

SbU.-O + £ M T > fi i (t) 

1-1 


T,(S,t) + £ M t) ^ Cs) 
1-1 


- * «£ /J r i(®> ,ln k i( T - 6 > i9 


^5 J* q(t»T)y i ({) dt (i - i, 2, . . .) 

pi 

r 0 / f,({»T)dt (i - o) 

J o 


ifs*. fi a <°) + y^to] (l - 1, a, • • 0 
(l + 5 ,.)y 0 e (i - 0 ) 


-J-*- f r f T f 1 5 (t,T)diCdT ) 2 

1 + n,- Jo JO Jo 


i:i:a i ^Th 


•s' Si iu ' T)(at>2 + - 1 + i sS iu ’ T)at 


- f ‘ q(E,T)dj + — i-^-U - l) f 1 5(t»T)at 

J 1 1 + »o J 0 


C jeosh jfk* cos + cob ^ cosh - 

h s U f[ ♦ .uh fc)U fi + «n ^ ~ — ^ ♦ "* Jp - ” 

~ Vi \ ’vV 1 ’ 1 \ .la fa + .ii* cos ^ + oo.h 


-Ch^- cosh ^ cos ^4 + cos ^ cosh ^4 + 

% £U * slab fc)(cos ^ ♦ cosh 

\slh ^ + slab ^ oo» +■ cosh J 


-Ckj^ 2 jcosh ^ sla ^4 + cos ^ slab |^4 + 

h (sin * slab i/ETVoos + cosh lfc{ * «”!» fat +. ala M ~ 

“ \ V V\ sla^iq+slah^I cos fc + cosh ^ } 



t - cos jh^4\ 
+ oosh y 


Frequency equations 


i Vk ainh ^k + Bin >lk co«h Vk + 5 q ^k(l + cos /k coah /k) - 0 
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TABIE I.- RESPONSE ® A UKEFOHM EIH4EHTARY BEAM TO A ffiNERAl LOAD - Continued 

(b) Cantilever beam 


Quantity 

Analytical expression 

y(ijT) 

00 

y s (s> T ) + ^ ^(t) y L (t) 

i=i 

A(5,t) 

00 

M s (£,r) + £ 0 ± (t) M 1 (£) 
i=l 

V(6,t) 

00 

V s (£,t) + £ ^(t) v ± (|) 
i-1 

^i(r) 

p (x) k. p T 

i -l- -1 / P-Ce) sin k.(r - 9)d0 

“i m i Jo 1 

Pi(-r) 

- 

^2 /q 1 5(£ ’ T) dS 

“i 

/j, 1 ^i 2( ^ dS ■ fc 

y e (s,T) 

r ! r' r s r‘ 5 (i,t)(.d , ‘ 

Jo Jo Ji Ji 

Ms(5,t) 

~ Si Si ^® ,T ^ d ^ 2 

V s (£,t) 

r £ _ 

- q.(£,T)d£ 

ytCs) 

/sin - slnh cosh fkli - cos ^s\ 

y sin fk^ +■ sinh cos ^k^ + cosh (/kj[ / 

Mi(£) 

(Sin + sinh ^5 cos /k^g + cosh fkit\ 

\ sin ^ + sinh fk^_ cos + cosh ^ j 

V ± (£) 

„ lr 3/2j cos ^is + oosl1 , sin - Binh i/^lA 

\ Bin fk± + sinh. cos \/k£ + cosh fk^ J 


Frequency- equation: 


1 + cos fk. cosh fk = 0 




















T6BIB I— RESPONSE CF A XB3FQRM ELEHKHTAEy BEAM TO A GENERAL LOAD - Concluded 


VJi 

ro 


(c) Simply supported Tjeam 


Quantity 

Analytical expression 

y<l,T) 

09 

1-1 

MU,t) 

CO 

flgU/r) + £ 1 Z x\(t) sin 1*6 
1-1 

V(£,t) 

oo 

? 0 (£,t) + ^ iV0j_(t) cos i*£ 
1-1 


-^(t) + 21 2 * 2 [" P 1 (0) sin 1 2 j£ 2 (t - 0)d0 

Jo 

P t (T) 

5(£,t) sin 1*£ d£ 

^ T ) 

r* r‘ r 5 r‘ 5( t ,w - 1 r 1 r s r s r s a.,w ♦ ^ - & r 1 r* 

Jo Jq Jo Jo J o Jo Jo Jo 0 Jo Jo 


- f f i(£j t )(4£) 2 + i f f 5(£,T)(d£) 2 
Jo Jo Jo Jo 

V b U,t) 

- [ t l(i,r)Q.i+ f f* q(t,T)(d£) 2 

Jo Jo Jo 


t I 


+T+TS+T 1 ® VOYK 





NACA TN 4244 


53 




Analytical expression 


y r (-r) + + £ jijto T^t) 










“.(S.-r) + £ * l ( t ) fijCf) 
1-1 


m 

v,(l, t) + £ ^(t) fi(t) 

i -1 


-54^4 + jjt s^e) .in ^(t - e)aa 


^5 5<{,t) f ± (t) at (i . i, 2 , . . .) 

r 1 

?o / f.(i,T)at (i - o) 

JQ 


(i * x ' £ ' • • 0 
(l + i^y 0 2 (1 - 0 ) 


r^/ 0 7o7c> T >*' (aT > 8 


7777 5Ct ' TK4{)k - <77 5({ ' T>(d,)E + rr^ly^ -‘)-(s-t + t)1 7 «»•*• 


777 * (£ ' T,(d£)5 - - t + i7 5<t ' T)d£ 


- 77 5(£ ' TKd£)E + - 8 + i) 7 5(£ ’ t)4{ 


- 5 (ijT) 4 i + _i_<i - i) i(i,T)d{ 


■ In S. »inll B» 

— p—i MMh Ojt 5 ^-^ co. - 

f s 2 - *A a ) r / ' 

4i__3_£-i^co.ti Oj + co. jco. Pjl - cortl Ojl + Pi 8 " {TJ ^ “l £ 
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TAB 12 II— RESPONSE CF A UNIFORM TIMOSHENKO BEAM TO A GENERAL LOAD - Continued 
(a) Symmetrical free -free bean vith. a concentrated aass - Concluded 
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TABLE H.- KESPOHEE OF A T3CTCEK IDDSEXOD BEAK ZO A QZKEBAL LOAD - Continued 
(b) Cantllem b east 


Qtmntlty 

Analytical expreaslcc. 

?(S,T) 

? B U,r) + £ fcOO y ± (|) 
1-1 

t<S,T) 

t*(t/0 + £ 
1-1 


MjCSjT) + Y W r) fi l Cl) 
1-1 

nt,r) 

f,(t,T) + £ am \cn 
1-1 

AW 

~%r^ + :r f p t c«y «iii i^Ct - e)iu 
*1 *1 J 0 

PjW 



//lAl) + ‘pi 2 -4 

f.UW 

r 1 r e r‘ r 1 «t,w - ^r* r 1 

Jo Jo Jl Jl Jo Jl 


r* r e r E 5(*, T )(a*j 3 

Jo Ji Ji 

^(e.r) 

'fill w, ' T>C4t)2 


-J* 5(|»t)4| 

fjCl) 

cjeoa p^t - coah. ajfc + Aj^aic. Pjft - ^ Jj* alnb Ojljjj 

ACt) 

o ^ ^i(°°* A* " ao * h A 6 ) ■ (“ ln At + A • 141,1 “it)] 

AW 

c (a 2 ■ A^» s )j^i( ,In A 1 + js£ ,Illh “it) ♦ («=» tjl + 7 t oorti *)] 


0 5^a(°°“ f±t * n 00,51 A 1 } - ( ,In A* - “it)J 


Pi 

sin Pi. - alnh. a± 

co* Pi + ^ coali 

mi 

“i 2 + A 2 *! 2 

A 2 - tiV 

WKM 

ki ji[-( l£ * a + 1 *4 + $S Z ' ‘to 2 ) + ^5 


HHI 




Frequency equation: 

2 + - ^jaln £ ainli a. + ^ + ^co* P cosh a - 0 
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TABLE II— RESPONSE CF A OBirOEM TDCSEGENKO BEAM TO A OBNffiAL LOAD - Concluded 
(o) Simply Buppcrted toean 






























Figure 1.- Positive distortions and positive internal forces and mome nts 
associated with a typical beam element. 
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Figure 2.- Grid scheme for traveling -wave numerical procedure . 










Figure 4.- Shear response of a uniform free-tree beam to a step load concentrated 

at § a 0. 
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Figure 5 .- Shear response 



of a uniform free -free heara. to a ramp-platform load 


concentrated at £ »= 0. 


O'! 

H 


NACA TN KZkh 



concentrated at | = 0. 
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Figure 7.- Shear response of a uniform free-free beam bo a triangular pulse load 

^nnofm+vr-n.+^rl at. P = 0. 

»- 3 


t+rSt HI V 0 VN 




o .2 A .6 .8 1.0 1.2 


T 

Figure 8,- Shear response V^ 7 t^ °f a imifcmii free-free heam to a blaBt pulse load 

concentrated at £ - 0. 







